Trap-imbalanced fermion mixtures 
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We analyze the ground state phases of two-component [a = {|, I}) population- and mass-balanced 
(iV| = TVj and m| = mj) but trap-imbalanced (ujy 7^ uy) fermion mixtures as a function of inter- 
action strength from the weak attraction Bardeen-Cooper-Schrieffer (BCS) to the strong attraction 
Bose-Einstein condensation (BEC) limit. In the BCS limit, we find that the unpolarized superfluid 
(UPS) fermions exist away from the central core of the trapping potentials, and are surrounded by 
partially polarized normal (PcrPN) fermions. As the interactions increase towards unitarity, we find 
that the central PcrPN core first transitions to a UPS, and then expands towards the edges until 
the entire mixture becomes a UPS in the BEC limit. 
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Ultracold atomic physics experiments with two- 
component fermion mixtures have enabled the study of 
novel superfluid and insulating phases which have not 
been possible in other systems. For instance, the tun- 
ing of attractive fermion-fermion interactions have per- 
mitted the ground state of the system to evolve from 
a weak fermion attraction Bardeen-Cooper-Schrieffer 
(BCS) limit of loosely bound and largely overlapping 
Cooper pairs to a strong fermion attraction limit of 
tightly bound and small bosonic molecules which undergo 
Bose-Einstein condensation (BEC) 0, 0, i, 0, H, @ . For 
mass- and population-balanced mixtures, in agreement 
with the early theoretical predictions 0, these experi- 
ments have shown that the BCS-BEC evolution is just a 
crossover. 

Recently, the ground state phase diagram of mass- 
balanced but population-imbalanced fermion mixtures 
have been theoretically analyzed showing that the 
BCS-BEC evolution is not a crossover but quantum 
phase transitions occur between normal and superfluid 
phases [1, 0|. In addition, phase separation between su- 
perfluid (paired) and normal (excess) fermions has been 
shown. Motivated by these predictions, there have been 
several experiments with mass-balanced but population- 
imbalanced fermion mixtures 110 . 11 1 . leading to an in- 
tensive theoretical activity [H, HI EI El, El EI El ■ 



Since exotic superfluid phases (i.e ., Fulde-Ferrell and 
Larkin-Ovchinnikov (FFLO) [H 21 1) can be potentially 
realized, imbalanced fermion mixtures are currently in- 
terest to many communities ranging from atomic and 
molecular to condensed- and nuclear-matter physics. For 
instance, the ground state phase diagram of mass- and 
population-imbalanced fermion mixtures have been re- 
cently analyzed showing quantum and topological phase 
transitions [II III, [II 
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In this manuscript, we analyze the ground state phases 
of two-component (cr = {T,l}) population- and mass- 
balanced (iVf = N± and m| = mjj but trap-imbalanced 
(w-f 7^ u>i) fermion mixtures as a function of interac- 
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FIG. 1: (Color online) Schematic diagrams showing shell- 
structures of two-component (cr = {|, J.}) population- and 
mass-balanced (Ny = Ny and ni| = my ) but trap-imbalanced 
(uiy > ujy) fermion mixtures for (a) non-interacting, and (b) 
weak, (c) intermediate and (d) strong attraction regimes. 
Here, uia is the trapping frequency of a-fermions. The colored 
regions correspond to unpolarized superfluid (UPS, red), un- 
polarized normal (UPN, white), partially |-polarized normal 
(P|PN, blue), partially ^-polarized normal (P|PN, green), 
and fully .[-polarized normal (FJ.PN, black) phases. 



tion strength, where ui a is the trapping frequency of a- 
component. The ground state involves very rich shell- 
structures consisting of unpolarized superfluid (UPS) 
and unpolarized normal (UPN) as well as partially a- 
polarized normal (PcrPN) and fully cr-polarized normal 
(FcrPN) fermions. Our results are schematically shown 
in Fig. [TJ and are as follows. In the BCS limit shown in 
Fig. DIb), we find that the UPS fermions exist only away 
from the central core of the trapping potentials, and are 
surrounded by PcrPN fermions. As the interactions in- 
crease towards unitarity, we find that the central PcrPN 
core first transitions to a UPS as shown in Fig. [TJc), and 
then expands towards the edges until the entire mixture 
becomes a UPS in the BEC limit as shown in Fig. [ljd). 

Pairing Hamiltonian: To obtain these results, we start 
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with the Hamiltonian density (in units of Ti = ks = 1), 
H(r) = ^W^WV.W - ff *f a (r)* T)i (r), (1) 



where V4-( r ) creates a pseudo-spin-cr fermion at posi- 
tion r, and tpj ^(r) = i/»|(r)^j(r) is the pair creation 
operator. In Eq. (JTJ), g > is the strengh of the 
attractive fermion- fermion interactions, and we defined 
K a {v) = -V 2 /(2m CT ) - /i CT (r) where ^ CT (r) = - U CT (r) 
is the local chemical potential. The global chemical po- 
tentials /v fixes the density n a — N a /V of each type 
of fermion independently, where N a is the number of 
cr-fermions and V is the volume. The term V a (r) = 
m a {ui 2 x x 2 + uj 2 y y 2 + lo 2 z z 2 )/2 corresponds to the trap- 
ping potential, which is assumed to be harmonic in space. 

In the momentum space, within the local-density (LD) 
approximation, the local mean-field (MF) Hamiltonian 
can be written as Hq(t) = S k)CT £k, CT (r)a ko .a k)(T ~ 

Aq(r) Ek( a I+Q/ 2 , T a -k+Q/2,i + h - c -"> + A q( r )/5> wher e 
Q is the center-of-mass momentum of individual Cooper 
pairs, and £k,<r(r) = e k)(T -/io-(r) with e ki<T = |k| 2 /(2m .). 
Here, Acj(r) is the local MF order parameter which is as- 
sumed to be real without loss of generality, and defined 
by A Q (r) = .9Ek( a -k+Q/2aflk+Q/2,T)' where (.) implies 
a thermal average. 

Self-Consistency Equations: The local MF Hamilto- 
nian can now be solved by using standard techniques [H, 
18l . l22j . The order parameter Aq (r) is determined by 



E 



f 1 l-/[ff k ,Q,T(r)]-/[ffk,Q,i(r)] 
\2e k 2£ k ,q(r) 



MV 

k " 

(2) 

where £k = ( e k.| + £ k,i)/2 is the average kinetic en- 
ergy, f(x) = l/[exp(x/T) + 1] is the Fermi function, 
and £ k , Q (r) = [£ iQ (r) + A^r)] 1 ^ with &lQ (r) = 
Kk+Q/2,T( r ) + C-k+Q/2,;(r)]/2. Here, £ k ,Q,<r(r) = 
^k,q(r) +s CT [^ k+Q/2:T (r) - £-k+Q/2,|(r)]/2 is the quasi- 
particle energy when s-f = 1 or the negative of the quasi- 
hole energy when Sj = — 1. Notice that, we ellimi- 
nate g in favor of the fcrmion-fcrmion scattering length 
ap via the usual regularization 1/g = —MV/i&nap) + 
E k l/(2ek), where M = 2m-]-m^/(mi + mjj is twice 
the reduced mass of f- and j-fermions. Eq. ([2]) has to 
be solved self-consistently with the number equations 
N a = J drn a (r) , where 



+ < Q (r)/[-Sk,Q,- CT (r)]}, 
is the local density of cr-fermions. 



(3) 
Here, 



*k,Q 



(r) = [1 + &,q(r)/£ k , Q (r)]/2 and i£ Q (r) 



[1 ~Ck.Q( r )/£'k,Q( r )]/2. Since the pseudo-spin symmetry 
is broken in population-, mass- and/or trap- unbalanced 
Fermi gases, one needs to solve all three equations 



self-consistently for all Q, and determine the value of Q 
which minimizes the free energy 0, 120I 21 1 . 

Having established the theoretical formalism, next 
we analyze the ground state phases of trap-imbalanced 
fermion mixtures as a function of of- For this purpose, 
we first discuss the numerical results and then provide 
analytical insight into the problem. 

Numerical Calculations: In this manuscript, we as- 
sume that the trapping potentials are isotropic in space 
such that Vo-(r) = m a Lo 2 r 2 /2 where r = |r|. In addi- 
tion, we do not explicitly consider the FFLO-like (Q ^ 0) 
superfluid phase 2^, 21 1, and limit the numerical calcu- 
lations to the Q — phase. However, this phase may 
be present in the weakly-attracting imbalanced fermion 
mixtures, but only in a narrow parameter space, as dis- 
cussed below. We also set N- 
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and 



consider two cases (a) u>i = l.lwj, and (b) u>-\ = 2o»j,. 
The numerical calculation involves self-consistent solu- 
tions of Eqs. {2} and (|3|) for Ao(r), n„(r) and /j, a . 
For instance, in Fig. [2j we show the local polarization 
6n(r) — u-f(r) — n^ir) as a function of of, characterizing 
the non-interacting, and weak, intermediate and strong 
attraction regimes. 




FIG. 2: (Color online) Local polarization 8n(r) — n-f(r) — 
nj, (r) [in units of fcf^|/(27r) 3 ] versus radius r (in units of Rf,i) 
is shown for population- and mass-balanced (JV-f = N± and 
77i| = mj,) but trap-imbalanced fermion mixtures, where the 
trapping frequencies are (a) wj = l.l^j. and (b) ui-\ = 2ujy. 
While the unpolarized regions are UPS, the polarized regions 
include both FFLO superfluid, and PaPN and FcrPN phases. 

We numerically find that the ground state involves very 
rich shell-structures consisting of UPS and UPN as well 
as PctPN and FcrPN phases, depending on the particu- 
lar value of ap as shown in Fig. [TJ To understand these 
shell-structures, next we analyze the non-interacting, and 
weakly- and strongly-attracting limits, which are analyt- 
ically tractable. 

Non-Interacting Fermion Mixtures: To understand the 
interacting trap-imbalanced fermion mixtures, it is use- 
ful to analyze first the non-interacting case when g = 
or ap — > — . In this limit, the mixture is in normal 
phase such that the superfluid order parameter vanishes 
at all space Ag(r) = 0, and that the global chem- 
ical potentials are identical to the global Fermi ener- 
gies Ha — ep^ a at zero temperature (T = 0). Thus, 
Eq. © reduces to n a (r) = (1/V)£ k /[&,„(»•)], and 
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at T = is given by n a (r) = (l/tO£ fc<fcj>(r) 1, 
where kp^(r) is the local Fermi momentum defined by 
e *> = k 2 Fa /(2m a ) = k 2 Fa {r) / {2 m<y ) + V a (r). This leads 
to n a (r) = mluil(R F (J - r 2 ) 3 / 2 / (6ir 2 ), where Rf,<j is the 
Thomas- Fermi radius of a- fermions defined by tp.cr — 
m a -uj 2 Rp a ./2, such that kF.a — m a uj a RF, a is the global 
Fermi momentum. Then, the number of cr-fermions is 
found by integrating n a (r) over r where r < Rf,<j, 
leading to N a = k Fa R F j48. Setting 7V T = N t gives 

Rf,]/Rf,i — [ m i^i/ 1 (m^u^)] 1 / 2 , and therefore, a trap- 
imbalanced fermion mixture can be realized when the 
condition m^ui^ ^ fn^i is satisfied. 
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FIG. 3: (Color online) Schematic (a) trap (V CT ) and (b) den- 
sity (n CT ) profiles are shown for non-interacting [1/ '(kp^ap) = 
— oo] population- and mass-balanced (TVf = iVj. and raf = 
m|) but trap-imbalanced (w-f > ui) fermion mixtures. 

When mfW-f > m^i, since i?i?T > Rf,i, there are 
more f-fermions near the center of the trap while J.- 
fermions are in excess near the edges as shown in Fig. |3] 
In addition, the local density of f- and J.-fermions are 
equal only at radius r c — RF,i[miuji(ui-[ — to \) / '(m^ui 2 — 
miLd 2 )] 1 / 2 , satisfying Sn(r c ) = 0. Therefore, in this case, 
the ground state corresponds to a PfPN for < r < r c , 
to a UPN for r = r c , to a P|PN for r c < r < R Fjh 
and to a FjPN for Rf,i < r < Rf,i- For instance, when 
= 2ojj, we find r c pa 0.58Rf,i and i?F,| ~ 0.71Rf,i, 
and these three phases can be seen in Fig. [2](b) when 
l/(fcF,i a f) = -co. The shell-structure of this case is 
schematically shown in Fig.QJa). 

Having finite and attractive fermion-fermion interac- 
tions changes this simple non-interacting picture dramat- 
ically, which is discussed next. 

Weakly- Attracting Fermion Mixtures: When g > 
or Of < 0, the normal mixture may become unstable 
against formation of Cooper pairs at some regions of 
the trap such that Ag(r) ^ 0. According to the BCS 
theory of superconductivity, this is always the case for 
chemical potentially balanced mixtures, 5/i(r) — [/ij (r) — 
/i|(r)]/2 = for all r, no matter how weak the g is. This 
suggests that, for an arbitrarilly small g, trap-imbalanced 
fermion mixtures first become unstable against superflu- 
idity at radius r c , where Sn(r c ) = and 5fi(r c ) = 0. 

In the weakly-attracting limit when g <C {ef.t? e fu}j 
the local order parameter is obtained from Eq. @, and 
is given by A (r) = (8/e 2 )/z(r) exp[7r/(2fc f (r)a F )], which 
is valid when /i(r) 3> Ao(r). Here, /i(r) = [pLf(r) + 
Hl(r)]/2 is the effective local Fermi energy and kp(r) — 



[2A//i(r)] 1 / 2 is the effective local Fermi momentum. No- 
tice that, Ao(r) has a maximum at the center of the trap 
and it vanishes towards the edges. For mass-balanced 
(rn| = mjj mixtures, a UPS phase with Q = minimizes 
the free energy when < \8^,{r)\ < 0.71Ao(r), however, 
further increase in \Sfi(r) | causes a first order transition to 
a PcrPN phase when \Sfj,(r)\ > 0.7lA (r) [H,^. There- 
fore, we find in the BCS limit that the UPS fermions are 
surrounded by the PcrPN fermions, and that they exist 
only away from the central core of the trapping poten- 
tials. For instance, these phases can be seen in Fig. Hfb) 
when l/(fci?,j,a_F) = — 1 or —0.7, and are schematically 
shown in Fig. QJb). Notice that, similar shell-structures 
have been also reported for purely mass-imbalanced mix- 
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tures 

As the interactions increase, Ao(r) increases at a faster 
rate near r = due to the faster increase in local fermion 
densities, which causes an additional first order transition 
from PcrPN to UPS at r = 0. For instance, when = 
2u>i, this transition occurs at 1/ '(kF^ap) ~ —0.65 as can 
be seen in Fig.[5{b), and the shell-structure of this case is 
schematically shown in Fig. [TJc). Notice that, this shell- 
structure does not occur with purely mass-imbalanced 
mixtures [24L 12a . l27j | . Further increasing the interactions 
towards unitarity, we find that the central UPS region 
expands towards the edges. 

In passing to the strongly- attracting limit, we make 
two comments. First, it is known that an FFLO- 
like superfluid phase with Ag(r) = A(r) exp[iQ(r)r] 
and Q{r) ~ 2AM\8n(r)\/kp(r) may also exist in a 
small parameter space when 0.71Ao(r) < \8^,(r)\ < 
0.75A (r) [13, EH . This phase resides between the UPS 
and PcrPN phases, and is separated from the UPS phase 
by a first order and from the PcrPN phase by a second or- 
der transition. Notice that, FFLO shells are not shown in 
Figs.QJb) andQJc). Second, the inclusion of fluctuations 
beyond the MF would reduce Aq(t'), and therefore the 
transitions discussed above are likely to occur at higher 
values of 1/ (kp.iap) than our MF predictions. While 
the weakly-attracting MF description is strictly valid for 
I/Qcf^clf) "C 0, it still serves as a qualitative estimator 
for the phase boundaries until l/(fcF,x a F) ^ —0.5. How- 
ever, this description can not be used for l/(A;F,j.a_F) ^ 0, 
which is discussed next. 

Strongly- Attracting Fermion Mixtures: In the strong 
fermion attraction (BEC) limit when g 3> {e_F,|, e_F,xl or 
o-f ^ 0, imbalanced fermion mixtures can be described by 
a mixture of weakly-repulsing Bose molecules and Fermi 
atoms, where the Bose molecules correspond to paired 
|- and j-fermions, and the Fermi atoms correspond to 
unpaired fermions 0,111. However, in population- 
balanced (iVj — Ni) mixtures, all \- and j-fermions are 
paired to form Bose molecules, and therefore the equa- 
tion of motion at T = is 

V 2 * B (r) 



2m i 



(4) 
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which is of the Gross-Pitaevskii form, where '^(r) = 
[M 2 a F /(87r)] 1/2 A (r) is the local BEC order parameter, 
Ms( r ) = 2/i(r) — £5 = (Xb— Vb(t) is the local chemical po- 
tential, Ubb = ^clbb/tub is the repulsive interaction, 
and tub = vti\ + is the mass of the molecular bosons. 
Here, [Lb = MT+Ml - e b is the chemical potential, Ve(r) = 
Vj-(r) + Vi(r) is the trapping potential, e& = — l/(Mop) 
is the binding energy, and clbb * af is the boson-boson 
scattering length of the molecules. Notice that, identifi- 
cation of Vb(t) = niB0J%r 2 jl as the molecular trapping 
potential leads to ojb = [(m^o; 2 + m|wj)/(m| + mj] 1 ' 2 , 
which is the effective trapping frequency felt by the 
molecular bosons. Therefore, we find in the BEC limit 
that the ground state of trap-imbalanced fermion mix- 
tures is the BEC of molecular bosons for the entire trap. 

These results are strictly valid for l/(fcf\i a F) ^S> 0, 
but they still serve as a qualitative estimator for the 
phase boundaries until 1/ {kp^a-F) > 1. For instance, 
when W| = 2a; j,, we find that the central UPS region 
expands, and the inner UPS and normal regions shrink 
towards the edges with respect to Fig.QJc), as l/(kF : yap) 
is increased. Notice that, the expansion of central UPS 
towards the edges is similar to the one observed with 
population-imbalanced fermion mixtures 0, How- 
ever, in our case, the normal regions vanish beyond a 
critical 1/ '(kp^ap), and the entire mixture becomes a 
UPS |30j. When toy = 2o>j, this occurs at l/(fc,F,i a F) ~ 
—0.1, and the shell-structure of this case is schematically 
shown in Fig. QJd). 

Having analyzed the ground state phases, next we dis- 
cuss briefly the validity of our results and also their ex- 
perimental realization in atomic systems. 

Experimental Realization: In this manuscript, we 
mainly rely on the LD, MF and isotropic-trap approx- 
imations. In LD approximation, the mixture is treated 
as locally homogenous, and this approximation is valid 
as long as the number of fermions is large [H, [3, , 
which is typically satisfied in atomic systems. In MF 
approximation, the superfluid order parameter is treated 
at the saddle-point level, and that the fluctuations are 
not included [22j. This description is qualitatively valid 
throughout the BCS-BEC evolution only at low tempera- 
tures [7[ , which can be reached in atomic systems. Lastly, 
in isotropic-trap approximation, the traps are assumed 
to be isotropic, while the atomic traps are typically elon- 
gated in one direction. However, the anisotropy of traps 
is not expected to affect the shell-structure of superfluid 
and normal phases other than causing shells to have el- 
liptical rather than circular cross-sections. 

In atomic systems, trap-imbalanced fermion mixtures 
can be realized in several ways. For instance, in the 
case of magnetically trapped systems, trapping two dif- 
ferent hyperfine states (f and j) of a particular atom (i.e., 
6 Li or 40 K) which have different magnetic moments (i.e., 
Ai] > Mi) corresponds to a situation where m-\ = mj 
and u>i > wy. Likewise, in optically trapped systems, 



asymmetrically detuning the laser frequency with respect 
to two hyperfine states may produce a state-dependent 
optical trap. Furthermore, trap-imbalanced fermion mix- 
tures can be naturally realized with two-species fermion 
mixtures [II [IS [III (i.e., 6 Li and 40 K) in both magneti- 
cally and optically trapped systems due to their different 
mass and also to hyperfine properties. 

Conclusions: We analyzed the ground state phases of 
two-component population- and mass-balanced but trap- 
imbalanced fermion mixtures as a function of fermion- 
fermion interactions. In the BCS limit, we found that the 
UPS fermions are surrounded by PctPN fermions, and 
exist only away from the central core of the trapping 
potentials. As the interactions increase towards unitarity, 
we found that the central PcrPN core first transitions to 
a UPS, and then expands towards the edges until the 
entire mixture becomes a UPS in the BEC limit. 

We thank C. A. R. Sa de Melo, P. S. Julienne, I. B. 
Spielman and R. Grimm for useful discussions. 
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